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Digital Filtering of Three-Dimensional Lower Extremity 
Kinematics: an Assessment 




i 



Errors in kinematic data are referred to as noise and are an undesirable portion of any waveform. Noise is 
typically removed using a low-pass filter which removes the high frequency components of the signal. The selection of 
an optimal frequency cut-off is very important when processing kinematic information and a number of techniques 
exists for the determination of an optimal frequency cut-off. Despite the importance of cut-off frequency to the efficacy of 
kinematic analyses there is currently a paucity of research examining the influence of different cut-off frequencies on the 
resultant 3-D kinematic waveforms and discrete parameters. Twenty participants ran at 4.0 m 9 s-l as lower extremity 
kinematics in the sagittal, coronal and transverse planes were measured using an eight camera motion analysis system. 
The data were filtered at a range of cut-off frequencies and the discrete kinematic parameters were examined using 
repeated measures ANOVA's. The similarity between the raw and filtered waveforms were examined using intra-class 
correlations. The results show that the cut-off frequency has a significant influence on the discrete kinematic measure 
across displacement and derivative information in all three planes of rotation. Furthermore, it was also revealed that as 
the cut-off frequency decreased the attenuation of the kinematic waveforms became more pronounced, particularly in the 
coronal and transverse planes at the second derivative. In conclusion, this investigation provides new information 
regarding the influence of digital filtering on lower extremity kinematics and re-emphasizes the importance of selecting 
the correct cut-off frequency . 
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Introduction 



Errors associated with the measurement 
of kinematic data can result from soft tissue 
artefact, improper digitization of retro-reflective 
markers and electrical interference (Winter et al., 
1974). These errors are typically referred to as 
noise and are an undesirable portion of any 
kinematic waveform. Noise is traditionally lower 
in amplitude and associated with a different 
frequency range than that of the true signal 
(Jackson, 1979) and can typically be removed 
using a low-pass filter; the objective of any 
filtering technique is not only to attenuate noise 
but also to leave the true signal unaffected 
(Winter, 1990). 



Several researchers (Enoka, 2008; 
Griffiths, 2006; Robertson et al., 2004; Winter, 
2005) introduced the concept of frequency content 
in retroreflective marker information. Filters are 
constructed so that different elements of the 
frequency spectrum are attenuated. When using a 
low-pass filter the cut-off frequency is selected so 
that the lower frequencies remain yet the higher 
frequencies are attenuated. The fourth-order zero- 
lag is frequently utilized in biomechanical 
analyses (Yu et al., 1999). Butterworth filters yield 
a weighted average of data points across the 
kinematic waveform. This method of averaging 
time points means that the filtered data lags 
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behind the raw data. To adjust for this time delay, 
the filter is passed through the data in the forward 
direction and then again in the opposite time 
order to yield data that is appropriately aligned 
temporally (Robertson and Dowling, 2003). In 
addition to attenuating any phase shift delay, the 
bi-directional filtering produces a sharper cut-off 
and is termed a fourth-order zero-lag shift (or 
zero-phase shift) filter (Smith, 1989). As the cut-off 
frequency increases, the influence of the filter on 
the data is reduced and the data will be similar to 
the raw signal, including some of the high 
frequency noise (Robertson and Dowling, 2003). 
Determining the most appropriate cut-off 
frequency is essential and necessitates knowledge 
of the signal's characteristics. 

The selection of a frequency cut-off is very 
important when filtering kinematic data. A 
number of algorithms exists which are designed 
to define an objective criterion for the 
determination of an appropriate cut-off frequency 
(Jackson, 1979). The first option is to perform at 
spectral analysis of the marker trajectories using a 
fast fourier transform (FFT) to examine the 
cumulative content of the signal in the frequency 
domain (Figure 1) (Giakas and Baltzopoulos, 
1997). Typically the choice of cut-off is taken as 
the frequency at which either 95 or 99% of the 
signal power is contained below (Sinclair et al., 
2012, 2013a, 2013b). 



A second option is to perform a residual 
analysis (Figure 2). Residual analyses examine the 
differences between the raw and filtered 
kinematic pattern over a pre-set range of cut-off 
frequencies (Wells and Winter, 1980). The term 
residual refers to the signal content that remains 
when then filtered data is subtracted from the raw 
data (Robertson, 2005). The geometric line A 
represents the best estimate of the noise residual 
and is positioned so that it follows the linear 
portion of the residual plot and intercepts the Y - 
axis at X (0 Hz). The decision regarding the cut-off 
requires a compromise between the extent of 
signal attenuation and the amount of noise 
allowed to pass through. Typically a horizontal 
line B is projected from X to intersect the residual 
plot at Y and the cut-off is selected at this 
frequency. 

Whilst it is recognised that the low-pass 
filter is the most commonly utilized smoothing 
procedure, the efficacy of the different methods 
for determining the optimal cut-off frequency 
remains unknown. The purpose of this 
investigation was to determine the effects of 
various low-pass filter cut-offs and optimal filter 
techniques on kinematic waveforms and discrete 
parameters. This study hypothesizes that the cut- 
off frequency will have a significant influence on 
the discrete 3-D kinematic parameters. 
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Figure 1 

Example FFT cumulative power analysis of marker information 
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Figure 2 

Example residual plot between filtered and unfiltered 
signals as a function of frequency cut-off 



Material and Methods 

Participants 

Twenty participants (sixteen males and 
four females) (age 26.33 ± 5.37 years, body height 
1.76 ± 0.12 m and body mass 75.5 + 8.6 kg) ran at 
4.0m-S" 1 ±5%. All were injury free at the time of 
data collection and provided written informed 
consent. Ethical approval for this project was 
obtained from the University of Central 
Lancashire School of Psychology ethics 
committee. 

Procedure 

Participants ran at 4.0 m-s -1 over a force 
plate (Kistler, Kistler Instruments Ltd., Alton, 
Hampshire) embedded in the floor (Altrosports 
6mm, Altro Ltd) of a 22 m biomechanics 
laboratory. Running velocity was quantified using 
Newtest 300 infrared timing gates (Newtest, Oy 
Koulukatu, Finland), a maximum deviation of 
±5% from the set velocity was allowed. Stance 
time was defined as the time over which 20 N or 
greater of vertical force was applied to the force 



platform (Sinclair et al., 2011). A successful trial 
was defined as one within the specified velocity 
range, where all tracking clusters were in view of 
the cameras, the foot made full contact with the 
force plate and no evidence of gait modifications 
due to the experimental conditions. 

Kinematic data was captured at 250 Hz 
via an eight camera motion analysis system 
(Qualisys Medical AB, Goteburg, Sweden). 
Calibration of the system was performed before 
each data collection session. Only calibrations 
which produced average residuals of less than 
0.85 mm for each camera for a 750.5 mm wand 
length and points above 3000 in all cameras were 
accepted prior to data collection. 

The marker set used for the study was 
based on the calibrated anatomical systems 
technique (CAST) (Cappozzo et al., 1995). In order 
to define the right; foot, shank and thigh retro- 
reflective markers were attached unilaterally to 
the calcaneus, 1st and 5th metatarsal heads, 
medial and lateral malleoli, medial and lateral 
epicondyle of the femur and greater trochanter. 
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To define the pelvis additional retro-reflective 
markers were placed on the anterior (ASIS) and 
posterior (PSIS) superior iliac spines. Rigid 
tracking clusters were positioned on the shank 
and thigh. Each rigid cluster comprised four 19 
mm diameter spherical reflective markers 
mounted to a thin sheath of lightweight carbon 
fibre with length to width ratios in accordance 
with Cappozzo et al. (1997). A static trial was 
conducted with the participant in the anatomical 
position in order for the positions of the 
anatomical markers to be referenced in relation to 
the tracking clusters, following which they were 
removed. 

Data processing 

Analysis of marker displacement data 
from the stance phase of the running cycle were 
examined using both residual and FFT analyses. 
Marker samples beyond toe-off were ignored, 
thus removing any non-stance data from the 
fourier reconstruction. Each subsample consisted 
of a total of 256 values, comprising N-adjusted 
marker values from a given contact phase and L = 
256-N zeros. The power (PW) of each component 
of the spectrum was calculated as the square of its 
amplitude (AM). Since the addition of zeros to the 
time-domain data reduces the calculated powers 
by a factor of N I (N+L), the inverse of this factor 
was applied to the power calculation to obtain 
true powers of the nonzero data. 

Thus, PW= (AM)\N+L) N- 1 

The residual analysis revealed an optimal 
cut-off frequency of 10Hz whereas the frequency 
analyses revealed frequencies of 25Hz 99% signal 
power and 10Hz 95% signal power. 

Hip, knee and ankle joint kinematics were 
quantified using Visual 3-D (C-Motion Inc, 
Germantown, USA) processed using raw data and 
also filtered using cut-off frequencies of 1Hz, 3Hz, 
5Hz, 7Hz, 10Hz, 15Hz, 20Hz, 25Hz using a zero- 
lag low pass Butterworth 4 th order filter. Lower 
extremity joint angles were created using an XYZ 
cardan sequence of rotations (Sinclair et al., 2012). 
All data were normalized to 100% of the stance 
phase then mean processed gait trial data was 
reported. Displacement specific 3-D kinematic 
measures from the hip, knee and ankle which 
were extracted for statistical analysis were 
1) ROM from footstrike to toe-off during stance 



and 2) peak angle during stance. First derivative 
angular velocity measures that were extracted 
were 1) peak stance phase angular velocity. 
Second derivative angular acceleration measures 
1) peak stance phase angular acceleration. 

Statistical analyses 

Descriptive statistics including means and 
standard deviations were calculated for each 
frequency cut-off. The statistical differences in 
each of the kinematic parameters as function of 
cut-off frequency were examined using repeated 
measures ANOVA's with significance accepted at 
the (p<0.05) level. Post-hoc pairwise comparisons 
were utilized for all significant main effects using 
a Bonferroni adjustment to control type I error. In 
addition to compare the similarity between the 
raw and filtered stance phase waveforms (i.e. 
extent of signal distortion/ attenuation) each of the 
mean filtered stance phase curves were contrasted 
against the normalized unfiltered waveform using 
intra-class correlations. 

Results 

Tables 1-3 present the intra-class 
correlation values between filtered and unfiltered 
waveforms. Figures 3-5 present the mean 
waveforms and discrete kinematic parameters 
obtained as a function of cut-off frequency. 

Angular parameters 

In the sagittal plane significant main 
effects were found for the magnitude of peak 
flexion p<0.01, if =0.54 and ROM p<0.01, if =0.89, 
post-hoc analysis showed that 1 Hz cut-off 
frequencies differed significantly from the 
unfiltered, 25, 20 and 15Hz conditions. In the 
coronal plane a significant main effect was 
observed for the magnitude of peak adduction 
p<0.01, r} 2 =0.69, post-hoc analysis showed that 1 
and 3 Hz cut-off frequencies differed significantly 
from the unfiltered, 25, 20 and 15Hz conditions. In 
the transverse plane a significant main effect was 
observed for the magnitude of peak external 
rotation p<0.01, r} 2 =0.74, post-hoc analysis 
showed that the 1 Hz cut-off frequency differed 
significantly from the others. 

In the sagittal plane significant main 
effects were found for the magnitude of peak 
flexion p<0.01, rf=0.80 and ROM p<0.01, r) 2 =0.53, 
post-hoc analysis showed that 1 and 3 Hz cut-off 
frequencies differed significantly from the others. 
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In the transverse plane a significant main effect 
was observed for the magnitude of peak internal 
rotation p<0.01, r] 2 =0.52, post-hoc analysis 
showed that 1 and 3 Hz cut-off frequencies 
differed significantly from the others. 

In the sagittal plane a significant main 
effect was observed for the magnitude of ROM 
p<0.01, r} 2 =0.53 post-hoc analysis showed that the 
1 Hz cut-off frequency differed significantly from 
the others. In the coronal plane a significant main 
effect was observed for the magnitudes of peak 
eversion p<0.01, n 2 =0.45, post-hoc analysis 
showed that the 1 and 3 Hz cut-off frequencies 
differed significantly from the others. 

Angular velocity parameters 

In the sagittal plane significant main 
effects were found for the magnitude of peak 
flexion velocity p<0.01, r] 2 =0.82 and peak 
extension velocity p<0.01, r] 2 =0.63, post-hoc 
analysis showed that 1 and 3 Hz cut-off 
frequencies differed significantly from the others. 
In the coronal plane significant main effects were 
found for the magnitude of peak abduction 
velocity p<0.01, r] 2 =0.62 and peak adduction 
velocity p<0.01, r) 2 =0.64, post-hoc analysis showed 
that 1 and 3 Hz cut-off frequencies differed 
significantly from the others. In the transverse 
plane significant main effects were found for the 
magnitude of peak internal rotation velocity 
p<0.01, r} 2 =0.80 and peak external rotation velocity 
p<0.01, r) 2 =0.72, post-hoc analysis showed that the 
1 Hz cut-off frequency differed significantly from 
the others. 

In the sagittal plane significant main 
effects were observed at the knee for the 
magnitude of peak flexion velocity p<0.01, r] 2 =0.56 
and peak extension velocity p<0.01, r] 2 =0.53, 
r] 2 =0.72, post-hoc analysis showed that the 1 Hz 
cut-off frequency differed significantly from the 
others. In the coronal plane significant main 
effects were found for the magnitude of peak 
abduction velocity p<0.01, r] 2 =0.63 and peak 
adduction velocity p<0.01, r} 2 =0.60 post-hoc 
analysis showed that 1 and 3 Hz cut-off 
frequencies differed significantly from the others. 
In the transverse plane a significant main effect 
was found for the magnitude of peak internal 
rotation velocity p<0.01, n 2 =0.58 post-hoc analysis 
showed that 1 and 3 Hz cut-off frequencies 
differed significantly from the others. 

In the sagittal plane significant main 



effects were observed at the ankle for the 
magnitudes of peak plantarflexion velocity 
p<0.01, r] 2 =0.61 and peak dorsiflexion velocity 
p<0.01, r} 2 =0.52, post-hoc analysis showed that 1 
and 3 Hz cut-off frequencies differed significantly 
from the others. Furthermore, significant main 
effects were found in the coronal plane for the 
magnitude of peak inversion velocity p<0.01, 
r] 2 =0. 77 and peak eversion velocity p<0.01, r] 2 =0.80, 
post-hoc analysis showed that 1 and 3 Hz cut-off 
frequencies differed significantly from the others. 
In the transverse plane significant main effects 
were found for the magnitude of peak internal 
rotation velocity p<0.01, r] 2 =0.62 and peak external 
rotation velocity p<0.01, r} 2 =0.64, post-hoc analysis 
showed that 1 and 3 Hz cut-off frequencies 
differed significantly from the others. 

Angular acceleration parameters 

In the sagittal plane a significant main 
effect was found at the hip for the magnitude of 
peak acceleration p<0.01, r} 2 =0.99, post-hoc 
analysis showed that each of the cut-off 
frequencies differed significantly from one 
another. In the coronal plane a significant main 
effect was observed for peak acceleration p<0.01, 
r} 2 =0.98, post-hoc analysis showed that each of the 
cut-off frequencies differed significantly from one 
another. In the transverse plane a significant main 
effect was observed for the magnitude of peak 
acceleration p<0.01, n 2 =0.71, post-hoc analysis 
showed that each of the cut-off frequencies 
differed significantly from one another. 

In the sagittal plane a significant main 
effect was observed at the knee for the magnitude 
of peak acceleration p<0.01, r) 2 =0.92, post-hoc 
analysis showed that each of the cut-off 
frequencies differed significantly from one 
another. In the transverse plane a significant main 
effect was observed for the magnitude of peak 
acceleration p<0.01, r) 2 =0.89, post-hoc analysis 
showed that each of the cut-off frequencies 
differed significantly from one another. 

In the sagittal plane a significant main 
effect was found at the ankle for the magnitude of 
peak acceleration p<0.01, n 2 =0.75, post-hoc 
analysis showed that each of the cut-off 
frequencies differed significantly from one 
another. In the coronal plane a significant main 
effect was observed for the magnitude of peak 
acceleration p<0.01, r] 2 =0.96, post-hoc analysis 
showed that each of the cut-off frequencies 
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differed significantly from one another. In the 
transverse plane a significant main effect was 
observed for the magnitude of peak acceleration 



p<0.01, r} 2 =0.97, post-hoc analysis showed that 
each of the cut-off frequencies differed 
significantly from one another. 
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Figure 3 



Angular angle displacement waveforms from the hip, 
knee and ankle as a function of cut-off frequency 



Table 1 



Intraclass correlations for the between filtered and unfiltered 
displacement waveforms as a function of cut-off frequency 
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Table 2 

Intraclass correlations for the between filtered and unfiltered 
first derivative waveforms as a function of cut-off frequency 
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Figure 5 

Angular second derivative waveforms from the hip, 
knee and ankle as a function of cut-off frequency 



Table 3 

Intraclass correlations for the between filtered and unfiltered 
second derivative waveforms as a function of cut-off frequency 
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Discussion 

The aim of the current investigation was 
to examine the effects of different filter cut-offs 
and optimal filter techniques on 3-D kinematic 
waveforms and discrete parameters. This study 
represents the first to examine the influence of 
different filters on 3-D displacement and derivate 
information. 

This study examined the influence of nine 
different cut-off frequencies on the displacement 
and also first and second derivative joint 
kinematic parameters. In support of the 
hypothesis the key observation from the current 
investigation is that the different cut-off 
frequencies have a significant influence of discrete 
kinematic parameters in all three planes of 
rotation across displacement and derivative 
waveforms. Qualitative examination of the 
resultant waveforms obtained as a function of cut- 
off frequency shows evidence of both under and 
over smoothing. This therefore emphasizes that 
notion that selection of an appropriate cut-off 
frequency is essential in order to achieve 
empirically meaningful findings and conclusions. 

Several techniques currently exist for the 
determination of an optimal cut-off frequency; the 
current investigation examined the influence of 
these techniques on the resultant frequency cut- 
off. Both the residual and fourier analysis at 95% 
of the signal power revealed cut-off frequencies of 
10 Hz. For joint kinematics from displacement 
and first derivative measures this cut-off appears 
to be stable as the waveforms were qualitatively 
free from noise yet comparative waveform 
analyses show that the signal is minimally 
attenuated. Therefore it seems that for the 
quantification of displacement and first derivative 
measures the residual and 95% signal power 
techniques are sufficiently reliable to provide a 
good estimation for the frequency cut-off. 
However, when using the spectral technique in 
which the cut-off frequency is taken at 99% of the 
signal power the same displacement and first 
derivative waveforms show clear evidence of 
noise. Therefore, whilst this cut-off frequency 
leaves the resultant waveforms minimally 
attenuated, it appears to be ineffective for the 
quantification of lower extremity kinematics even 
at the lower derivatives. 

When examining lower extremity 
kinematics of second derivative angular 



acceleration curves it is evident that the digital 
filtering can have a much more pronounced effect 
on the kinematic waveforms compared to those of 
the lower derivatives when the same cut-off 
frequencies are used. This suggests that the signal 
power of the higher derivative information is 
contained at a much higher frequency. Like the 
majority of kinematic analyses the optimal cut-off 
frequency was calculated by considering only the 
displacement data which was subjected to be 
spectral and residual analyses (DAmico and 
Ferrigno, 1990; Jackson, 1979; Winter et al., 1974; 
Yu, 1989). The higher derivatives were quantified 
by filtering the displacement data and then 
differentiating in the time domain (Damico and 
Ferrigno, 1990). However as mentioned above, the 
frequency characteristics of a signal are likely to 
be different for derivatives of a different order 
(Yu, 1989). Therefore, it appears that when 
calculating second derivative kinematics through 
successive differentiation of displacement data, 
there is amplification of the noise inherent in 
apparently smooth displacement curves. This 
suggests that both residual and spectral analyses 
of displacement data are not adequate for the 
determination of cut-off frequencies for higher 
derivative kinematics. 

Time-frequency analysis in which the 
signal is considered in both time and frequency 
domains simultaneously has been proposed as an 
alternative to traditional low pass filtering for 
derived data. Time frequency analysis techniques 
typically consist of short-time Fourier transforms 
and wavelet or Wigner distributions to localize 
the frequency content of the signals (Stankovic, 
2000). However issues do exist with time- 
frequency based filtering techniques. One of the 
failings of the short-time fourier transform is its 
fixed resolution in which the width of the 
windowing function relates to how the signal is 
represented. A longer window will give a good 
frequency resolution but result in a poor temporal 
resolution, whereas a narrow window function 
gives good time resolution but poorer frequency 
resolution. Furthermore, wavelet analyses can 
also be affected by pseudo-gibbs artefacts which 
may produce large oscillations in the time signal. 
Hatze (1981) presented a digital filtering 
technique in which the optimal filtering factor is 
quantified separately for each derivative domain. 
Therefore, instead of first filtering and then 
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differentiating the displacement data, the data are 
first differentiated and then filtered. This 
procedure assures a truly optimal regularisation 
for each of the derivatives (Hatze, 1981). Vaughan 
(1982) proposed that the displacement data be 
filtered using a cut-off frequency different from 
the optimum cut-off frequency derived from 
displacement data. Both Antonsson and Mann 
(1985) and Giakas and Baltzopoulos (1997) 
recommend separate filtering in each derivative 
domain following estimation of the optimal cut- 
off using spectral analysis although they did not 
indicate at what level of the frequency spectrum 
the cut-off should be implemented. The efficacy of 
these alternative differentiation and filtering 
techniques has yet to be evaluated although an 
understanding of which strategy is more effective 
will offer a means to improve the performance of 
existing and future automatic smoothing 
techniques based on digital filters. 

A further interesting observation of the 
current investigation is that regardless of cut-off 
frequency kinematic waveforms be more 
attenuated by the filter in the coronal and 
transverse planes particularly at the first and 
second derivatives. This suggests that the smaller 
more finite rotations outside the sagittal plane 
may be associated with a different frequency 
content in predominantly single planar motions 
such as running. Therefore, it may be that a 
different cut-off frequency is required for the 
different planes of rotation although additional 
work is necessary to confirm this. This 
observation may be less pronounced in less 
sagitally dominated movements such as walking 
or in movements in which the sagittal plane is not 
the prevailing rotation. 

A potential limitation of traditional 
fourier analyses in kinematic analyses is that the 
examined signal is truncated to just the stance 
phase of the gait cycle. Fourier analyses assume 
that the finite data represents one period of a 
periodic signal (Walker, 1996). For the FFT, both 
the time and frequency domain are circular 
topologies, so that the two endpoints of the time 
domain signal are interpreted as though they are 
connected together (Welch, 1970; Harris, 1978; 
Gautam et al., 1995). Therefore, the finiteness of 
the sampling record may result in a truncated 
waveform with different spectral characteristics 
from the original continuous-time signal. This 



truncation of the signal can introduce sharp 
transition changes into the measured data, a 
phenomenon known as spectral leakage. These 
sharp transitions are discontinuities and can 
influence the resultant frequency spectrum. To 
minimize this effect, we can apply a window 
function to the measured signal in the time 
domain. This will make the endpoints of the 
waveform meet and therefore result in a 
continuous waveform without sharp transitions. 
Windows are useful in reducing spectral leakage 
when using the FFT for spectral analysis, 
however, because windows are multiplied with 
the acquired time-domain signal, they introduce 
distortion effects of their own. The windows 
therefore change the overall amplitude of the 
signal. 

A further consideration of the current 
investigation is that a single cut-off frequency was 
applied to all of the markers. This is commonplace 
in kinematic analyses but does raise concerns as it 
is probable that more proximally positioned 
markers will be associated with a different 
frequency content than those placed more 
distally. A select number of researchers have used 
different frequency cut-offs for markers 
positioned at different locations (Mills et al., 2011) 
although this would drastically increase the 
duration required for data processing it may 
provide a more meaningful data set. Future work 
should consider contrasting the kinematic 
parameters and waveforms when an optimal cut- 
off has been determined for specific markers 
against a universal frequency cut-off. 

This study examined the influence of a 
low pass digital filter with different cut-offs on 
the same data set, the rationale behind this is that 
this smoothing mechanism is the most widely 
utilized in kinematic analyses (Winter et al., 1974; 
Yu et al., 1999). However, there are many different 
smoothing methods available in the literature in 
addition to conventional digital filtering; such as 
cubic splines (McLaughlin et al., 1977; Soudan and 
Dierkx, 1979; Zernicke et al., 1976), quantic splines 
(Wood and Jennings, 1979), polynominal 
smoothing (Pezzack et al., 1977) and fourier series 
(Andersson and Bloomfield, 1974). It appears that 
it would be of practical significance for further 
analyses to be conducted to examine the effects of 
different smoothing techniques at different cut-off 
frequencies on the same data set. This would 
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further enhance knowledge regarding the optimal 
filter mechanism/ cut-off frequency for the 
representation of movement analyses. 

In conclusion, the results of the current 
investigation suggest that different cut-off 
frequencies can have a significant influence on the 
discrete kinematic parameters. Inspection of the 
waveform analysis also suggests that for the 
quantification of 3-D kinematic parameters for 
displacements and from the lower derivatives that 



residual and fourier analyses using displacement 
information from lower extremity markers are 
able to satisfactorily produce an empirically 
meaningful cut-off frequency estimate. However 
future work is necessary in order to determine the 
most effective filtering technique when examining 
lower extremity kinematics at the higher 
derivatives. 
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